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Abstract 

We discuss the effect of an homogeneous flow in the aggregation process of 
colloidal magnetic particles at moderate concentration. Situations in which 
the presence of flow acts in favor of the chaining process: particles assemble 
into chains larger than the ones emerging in the absence of flow, under the only 
influence of an externally imposed field, have been analyzed. The results we 
obtain follow from the analysis of the pair correlation function which, owing 
to the potential character of the flow we consider, can be interpreted in terms 
of a Boltzmann-like stationary distribution function. To render the influence 
of the flow on the resulting structures explicit, we study the particular cases 
of axisymmetric and planar elongational flows. 
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I. INTRODUCTION 



One of the most remarkable features of the many-particle regime of suspension of mag- 
netic colloids is the formation of structures -mainly chains or branched structures- which 
may in general evolve over timeS. The reason for the occurrence of this phenomenon is the 
presence of dipolar interactions among particles. In the absence of an external field, the 
form of the emerging structures depends basically on the ratio between dipolar and ther- 
mal energies. When dipolar energy is dominant, alignment of particles is favored making 
the chain the most likely structure. The tendency to disorder, promoted by thermal mo- 
tion, may destroy that alignment resulting in the appearance of more isotropic aggregates. 
These structures are in fact fractal objects and have been classified in terms of the fractal 
dimension!'!. A fractal dimension close to one indicates the presence of chain-like aggregates 
whereas higher values of that quantity is the signature of DLA magnetic aggregates!. 

When an external magnetic field is applied, the dipoles are oriented along the direction 
of the field thus favoring the formation of chain-like aggregates. This strongly anisotropic 
situation becomes manifest through the form of the pair correlation, revealing that the 
underlying structure is the chami. 

The earlier stages of the chaining process are characterized by the existence of a scaling 
regime for the size of the chain-like aggregates. In the absence of flow, the scaling behavior 
for the size of the chain S ~ t a has been found with an exponent a = 1/2, obtained under 
the approximation of spherical cluster i. By considering the more realistic situation in which 
the diffusion process is mediated by anisotropic hydrodynamic interactions among particles, 
one obtains logarithmic corrections to the power laws providing the effective value of the 
exponent a e ff = 0.61 i, which agrees with experimental results I. At longer times, one 
has to consider the breaking of the chains when they become larger than a characteristic 
length. This feature is responsible for the existence of a crossover to a stationary regime 
characterized by the average size of the chains. It is in this regime in which we want to 
analyze the influence of the flow. 
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The situations reported previously occur when aggregation phenomena take place in a 
quiescent liquid. In many cases of interest, i.e. porous media, active dampers, rheological 
devices to control the viscosity, or in general in rheology of field-responsive systems, aggre- 
gation proceeds under conditions creating a stationary flow. The presence of these chain-like 
structures in turns, affects the rheology of the systemi _ 0. The question is to what extent 
the flow constitutes an important factor to determine the shape of the emerging structures. 
This is precisely the problem we address in this paper. We study aggregation of the mag- 
netic grains under an elongational flow. Owing to the potential nature of that flow, for 
which the stationary distribution function of a pair of particles is Boltzmann-likeffil, we can 
perform an analysis similar to the one carried out for aggregation in a quiescent liquid!. The 
knowledge of the pair correlation function or equivalently of the second virial coefficient will 
allow us to derive the equation of state in the presence of flow and to compute the mean 
size of a chain. 

We have organized the paper in the following way. In section 2 we describe the system and 
discuss the nature of the interactions between particles. Section 3 is devoted to calculate the 
pair correlation function and the related second virial coefficient corresponding to different 
situations concerning the relative orientation of the field and the external flow. The average 
size of a chain relative to the different situations discussed in the previous section is calculated 
in section 4. We extend the results of previous theories, expressing that quantity as a function 
of the magnetic energy between particles, to the case of an externally imposed elongational 
flow. Finally, in the last section we summarize the effects that the flow has on the kinetics 
of chaining processes. 

II. INTERACTING PARTICLES IN EXTERNAL FLOW 

We consider a colloidal suspension of magnetic dipoles in a non-polar liquid phase. The 
particles bear permanent magnetic moments rigidly attached. The one for the i-th particle 
is 
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mi = mRi , (1) 

where Rj is the unit vector along the direction of m; and m is the constant magnetic moment 
strength. The system is under the influence of both a constant external magnetic field H, 
high enough to induce alignment of the magnetic grains, and a stationary elongational flow, 
described by the velocity field 

Vi = AijXj , (2) 

where Aij are the components of the velocity gradient tensor, and Xj is the j-th component 
of the position vector. 

The flow induces on a pair of bounded particles a stretching, with associated energy 

Ufi = AfiMffi- (3) 

which depends on the components of the vector r joining the center of the j-th and i-th 
spheres. Here £ is the Stokes friction coefficient. 

The potential energy of a pair of particles splits up into the following contributions 

U = U fl + U mag (4) 

where apart from the potential associated to the flow defined in Eq. (|3]) we have introduced 
the energy U mag resulting from interactions between two magnetic dipoles a distance r apart 

If 3 1 

Umag = -3 j mi • m 2 - -(mi • r)(m 2 ■ r) I . (5) 

As follows from expressions (§)-(H), the energy of a pair of particles, and consequently the 
kinetics of the chaining process, depends on both the velocity gradient, and the orientation of 
the magnetic field with respect to the symmetry axis of the flow. To analyze that dependence 
explicitly, two situations will be under scrutiny. First, we will consider a planar elongational 
flow, which can be generated by four rolling rolls placed at the corners of a square, for which 
the nonzero components of the velocity gradient tensor are given by§ 
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A 12 = A 21 = A/2 , (6) 

and second an axisymmetric elongational flow, representing the flow through a pore, for 
which the nonzero components of the velocity gradient tensor are given by III 

An = A 22 = -A/2 ; A 33 = A , (7) 

with A/2 being the rate of elongation. 

III. THE PAIR CORRELATION FUNCTION AND THE SECOND VIRIAL 

COEFFICIENT 

For concentrations pertaining to the semi-dilute regime, the ferromagnetic particles tend 
to form chains along the direction of the imposed magnetic field. The system then reaches 
an equilibrium state in which the distribution function of the relative position of a pair of 
bounded particles is given by the corresponding Boltzmann distribution!. The flow acts as a 
perturbation such that the resulting stationary distribution is also a Boltzmann distribution 
in which the potential is given by Eq. (Q). It has been shown that this distribution is the 
stationary solution of the corresponding Fokker-Planck equation!. 

The pair correlation function of the suspension is given by@ 

g(r) = exp (—^=) (8) 



k B T, 

where T is the temperature and ks is the Boltzmann's constant. Its knowledge, allows us to 
derive the expression for the second virial coefficient, which is defined through the expression 

B{T) = -\ f (g(r)-l)dh. (9) 
2 Jn 

Here Q represents the accessible volume for a particle of the chain. 

In this section, the pair correlation and the second virial coefficient will be computed 
in the situations corresponding to the two types of flow introduced previously, taking into 
account different orientations of the magnetic field. For the planar elongational flow we will 
analyze the following cases: 



a) The magnetic field is oriented in the direction forming an angle 7r/4 with respect to the 
x-axis (see inset of Fig. 1). In this case, the pair correlation function (^j) can be written as 
follows 



9(r) 



exp {fx 2 sin 2 9 sin(2y?) + Ax~ 3 § sin 2 9(1 + sin(2y?)) - 1 } ; x > 1 



(10) 







x < 1 



where x = r/d, with d being the diameter of a dipole, and < 9 < n and < (p < 2n 
are the polar spherical angles. Moreover, we have introduced the dimensionless parameters 



A 



nr 



''/d^ksT, which corresponds to the ratio between the dipolar and thermal energies, 



and the Peclet number, P = ^k^T ' ex P ressrn g the ratio between hydrodynamic and thermal 
energies. Through this paper we will assume that A > P. To arrive at this formula, we have 
taken into account the hard core interaction between two dipoles. 

By substitution of Eq. ( |TD| ) into Eq. (Q) it follows the asymptotic expression of B(T) for 
large values of A 

-vrd 3 exp(2A + P/2) 



B 



(11) 



(6A-P)^/(3A + P)(3A + P/2) 

In obtaining Eq. fllT]) we have taken into account that the main contribution to the integral 
in Eq. (|9|) corresponds to r ~ d, and to the region around the directions 9 pa tt/2, and 
ip pa it/ '4 or ip pa 5ir/4. Notice that the dipoles link up preferably in the direction of the 
magnetic field, where the magnetic energy is attractive. Details about the derivation of Eq. 
(pT|) have been given in the Appendix. 

b) The magnetic field is now oriented parallel to the x-axis (see inset in Fig. 2). The pair 
correlation function is given by 



9{r) 



exp {f x 2 sin 2 9 sin(2<£>) + Ax 3 (3 sin 2 9 cos 2 (p — 1)| ;x>l 



(12) 



x < 1 



In this case, the asymptotic expression of B(T) is 



B 



71 



d* 



18A S 



exp(2A + P 2 /12A) ( 1 



V 3A 



(13) 



This result follows after consideration of the main contributions to the integral in Eq. @, 
which now correspond to r « d, 9 Pd tt/2, and (p « or y9 « 7r. 

c) The magnetic field is perpendicular to the plane of the flow. In this case, for the pair 
correlation function (BD we obtain 



9{r) 



exp {fx 2 sin 2 6sm(2ip) + Ax 3 (3 cos 2 cos 2 y? — 1)| ; x > 1 



(14) 



; x < 1 

Analogously, the asymptotic expression of -B(T) is now given by 



B 



nd 3 exp (2 A) 
6A J9X 2 - P 2 /4 



i + ^i n-i A 



3A 



2 9A 2 -P 2 /4 y 



(15) 



In deriving this expression we have used the approximation that the integral in Eq. @ 
extends only over the regions around 9 « 0, or 9 « 7r, and for r « d. 

For the case of an axisymmetric elongational flow defined through Eq. (0) we will consider 
that the magnetic field is oriented parallel to the z-axis. In this case, the pair correlation 
function (13) is given by 



g(r) 



exp { (f x 2 + Ax~ 3 ) (3cos 2 0- 1)} ;x>l 



(16) 



; x < 1 

By substitution of Eq. (|16D into Eq. (f|), we can obtain the asymptotic expression of B(T) 

exp (2A + P) 



B 



(17) 



(6A-2P)(3A + |P) ' 

where the integration domain of the corresponding integral reduces to the regions around 
the directions 9 ~ or 9 w 7T, and r ~ d. 
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IV. THE AVERAGE SIZE OF A CHAIN 



Our purpose in this section is to ascertain how the size of a chain is modified by the 
imposed flow. The average size of a chain composed of magnetic grains was first calculated 
in Ref. i as a function of the interaction parameter A. That approach consists in assuming 
that the chains of different sizes constitute a mixture of ideal gases in which the different 
components are identified with the chains having the same number of units. The pressure 
is then given through its expression for a mixture of ideal gases 



CO 



p = k B TY,Pk, (18) 
fe=i 

with pfc being the concentration of chains containing k dipoles. One defines the average size 
N of a chain as 

where p = SfcLi kpt is the initial concentration of dipoles. Using in this expression the virial 
expansion for p up to second order in p 

P 



p(l + pB), (20) 



the average size is given by 



N= r-r^T, ( 21 ) 

1 - p\B{T)\ ' V ; 



where we have taken into account the negative character of the second virial coefficient. This 
approach only provides the correct size of a chain for small values of A. At larger values 
of that parameter N diverges thus indicating that the actual kinetics of the process is not 
satisfactorily described in that domain. In order to avoid this imbalance between theory 
and experiment, a different approach, based essentially on the existence of a dynamical 
equilibrium among different species, was proposed 0. The kinetic process for which chains 
are formed is modeled as a set of coupled binary chemical reactions, in which two chains 
composed of k — n and n particles may 'react' resulting in a chain of size k or vice versa. 
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After defining the chemical potential of each species or energy per chain, and assuming the 
condition of chemical equilibrium one achieves a recursion relationship for p^ in terms of pi, 
which allows one to express the average size as 

N= , ^ . (22) 
'l + 4p\B\ - 1 



where p = n( p^ 6 ,with <p being the initial volume fraction of particles. Both theories report 
an increase of N as a function of the interaction parameter. The main difference lies in the 
fact that the second theory predicts the existence of finite chains for values of A larger than 
the ones corresponding to the first theory, for which N diverges. Divergence of the length 
of the chain occurs at A ~ 2 in Eq. (|2l|) , and at A ~ 6 for the case of Eq. (p2|) . Thus, 
to analyze the behavior of iV as a function of the elongational rate we will use the second 
approach, exhibiting a wider domain of physical values of A. 

The expressions of the virial coefficient B obtained for the different situations examined 
in section 3 can now be used in Eq. (B3) to compute the corresponding sizes of the chain. 
In case a), the averaged size of a chain becomes 



N =m exp(2A + P/2) 



1 + 2 4^. exp(2A + P/2) 



(6A-P)\/(3A + P)(3A + P/2) [ \ (6A - Pty(3A + P)(3A + P/2) 

(23) 

This result has been depicted in Fig. 1, from which one observes an enhancement in the 
length of the chain with respect to the case at zero elongation. 
In case b), N is found to be 

N = ^ ex P (2A + P 2 /2A) j-1 + ^1 + i£ exp(2A + P*/2\) J . (24) 

and has been represented in Fig. 2. This figure clearly shows corrections smaller than the 
ones observed in the previous case. This fact may be interpreted as a consequence of the 
loss of symmetry of the system inherent to the orientation of the magnetic field with respect 
to the flow. From the inset in this figure, we can realize that the effect of the flow is to drag 



away the particles towards directions different to that of the magnetic field. For this reason 
a stronger magnetic field, able to balance out the hydrodynamic forces, must be applied, 
and consequently the effect of the flow is less important than in the previous case. For case 
c) in which the averaged size of the chain is given by 



0exp(2A) 



A^9A 2 - P 2 /4 



-1 + 



1+ 7 - > 

\ Ay9A 2 - P 2 /4 



20exp(2A) 



(25) 



the flow does not introduce any significant correction to the result at zero elongational rate. 
This fact is a consequence of the particular orientation of the magnetic field, perpendicular 
to the plane of the flow. Due to the high strength of the magnetic field, the chains are 
oriented in a direction parallel to this field, and the chainig process results unaffected by the 

flow, irrespective of its nature. Finally, when the flow is axisymmetric we obtain 

-l 



1 + 24 J, - "Pff + f -V . (26) 



(6A - 2P) (3P/2 + 3A) 1 \| r (6A - 2P) (3P/2 + 3A) 

This result has been plotted in Fig. 3 and corresponds to the situation in which the presence 
of the flow influences to a greatest extent the chaining process. Indeed, for the same condi- 
tions as in case a), the correction to the value of N found in the absence of flow practically 
doubles the one obtained in that case. Now, due to the combined action of both magnetic 
field and flow, the grains become aligned along the symmetry axis of the system, thus favor- 
ing the chaining process. By comparison of the inset of this figure with the corresponding 
one in Figs. 1 and 2 one can realize that this is the more symmetric situation. 

From the cases analyzed previously one infers a common behavior: the presence of the 
flow results in an enhancement of the length of the chain when the elongational rate increases. 
In the case of the axisymmetric flow, when both magnetic field and elongation axis are 
aligned, the corresponding effects on the grains act in favor of the formation of a chain. For 
this reason, is in this situation in which larger chains may emerge. 
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V. CONCLUSIONS 



In this paper we have analyzed the influence that the presence of an external flow exerts 
on the distribution of sizes of chain-like aggregates in magnetic colloids. Due to the poten- 
tial nature of the flow we have considered, the potential energy of the resulting structures 
contains a contribution coming from the flow in addition to the dipolar magnetic energy. 
This fact has implications in the form of the pair correlation function and consequently in 
the value of the second virial coefficient. To assess that influence we have implemented typi- 
cal situations in which the magnetic colloid is subjected to conditions inducing elongational 
flows. The mean size of the chain can be calculated after computation of the second virial 
coefficient. To this purpose, we may apply different theoriesHlil originally proposed for a 
quiescent liquid. We have adopted the theory introduced in Ref. providing results valid 
for larger values of the A parameter. 

The potential character of the flow leads to a Boltzmann-like stationary distribution 
function for the configuration of bounded pairs of particlesEl, which allows us to perform 
an equilibrium-like analysis of the aggregation of particles in the presence of flow. Thus, 
we can obtain the equation of state that depends on the elongational rate, the parameter 
characterizing the flow. An equation of state for a liquid in the presence of flow in which 
the pressure depends on the shear rate was also obtained from thermodynamic grounds and 
simulations in Ref. 0. 

Among the cases we have discussed, we note that the most favorable situations for 
chaining are those in which the magnetic field lies along the directions corresponding to the 
minima of the elongational potentials. In these cases, magnetic field and flow cooperate to 
create stable chains. 

An estimate of the dimensionless parameters appearing in our theory, useful for possible 
experimental implementations, can be carried out by assuming particles of mass density 1 
g/cm 3 having a radius ranging from 10~ 6 cm to cm, suspended in water with viscosity 
10 -2 poise. We have also considered the elongational rate ranging from 10 s _1 to 100 
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and the magnetization of the particles of 10 Oe to 100 Oe. In all cases, the temperature has 
been assumed to be 300 K. Under this conditions, we can estimate the intervals 5 x 10~ 3 < 
A < 6 x 10 2 , and 2 x 10~ 4 < P < 2. The value of the different quantities corresponding to 
Figs. 1-3 are such that the dimensionless parameters belong to those intervals. 

Our main conclusion is that the external flow may introduce significant changes in the 
kinetics of the process resulting in the formation of chains larger than the corresponding 
ones obtained in the absence of flow. Our analysis has then revealed that the external flow 
constitutes an important factor in the formation of chain-like structures in magnetic colloids. 
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APPENDIX 

In this appendix we indicate how to compute the virial coefficients obtained in section 
3. For instance in case a), we assume that the dipoles forming the dimer are one close to 
other in such a way that x approximates to x ~ 1 + £, with £ << 1. In this situation the 
more important contribution to the integral @ comes from to directions: #~7r/2, ~ 7r/4 
and 7r/2, 57r/4, in which the magnetic energy dominates. Expanding the expression for the 
factor $ = —U/kT in those directions one obtains 

$ ~ 2A + P/2 - (3A + P/2) (6 - vr/2) 2 - (3A + P){<f>- tt/4) 2 - (6A - P)£ (27) 

and 

$ ~ 2A + P/2 - (3A + P/2){6 - vr/2) 2 - (3A + P)(<f> - 5vr/4) 2 - (6A - P)£ , (28) 
respectively. 
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To compute the virial coefficient one rewrites expression 



I r d 6 
B ~ -- / ci 3 rexp($) ~ -— exp(2A + P/2)I 1 I 2 I 3 (29) 

where I\, I2 and I3 are the following integrals 



00 



h = I de(l + 20 exp[-(6A - P)£] = - J— (l + -1- ) (30) 







6A - P V 6A - P 



h = jTdfl sin 0exp[-(3A + P/2)(6 - tt/2) 2 ] ~ ^ — ^ exp[-l/4(3A + P/2)] (31) 



and 

/ 37r/4 /■ 7-71-/4 

#exp[-(3A + P)(0-7r/4) 2 ] + / #exp[-(3A + P)(<f> - 5tt/4) 2 ] ~ 2. 
-tt/4 J3tt/4 ■ 



7T 



-tt/4 J3ir/4 V 3A + P 

(32) 



The expression (pTTJ ) for the second virial coefficient comes straightforwardly by combining 
equations (p9|)-(|32j). The virial coefficients for the remaining cases can be computed in a 
similar way. 
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FIGURE CAPTIONS 

Figure 1. Averaged size of a chain as a function of A for different values of P, and for 
= 0.14. The direction of the magnetic field form an angle of 45 degrees with the x-axis. 
The inset shows the corresponding orientation of the chain. 

Figure 2. Averaged size of a chain as a function of A for different values of P, and for 
4> = 0.14. The magnetic field lies parallel to the x-axis. The inset indicates the orientation 
of the chain. 

Figure 3. Averaged size of a chain as a function of A for four different values of P, and 
for 4> = 0.14. Now the magnetic field is parallel to the axis of the elongational flow. The 
inset also shows the position of the chain. 
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